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We study Goldstone boson (or axion) radiation from global C/(l)-strings 
^ ■ making direct quantitative comparisons between numerical scalar field the- 

I ory simulations and linearized analytic calculations in the dual antisym- 

metric tensor representation. Concentrating on periodic long string con- 
figurations, we show excellent correspondence for both the amplitude and 
spectrum of radiation. We also find good agreement with a linearized an- 
i alytic model for radiation backreaction, demonstrating that damping due 

G^ ■ to the lowest n=2 harmonic is effective in strongly suppressing higher har- 

monics. This work establishes the validity of the Kalb-Ramond action in 
p i' describing low-energy global string dynamics, along with the associated 

self-field renormalization of the equations of motion. We also describe the 
different nature of massive radiation from high curvature regions and the 
^ I collapse and annihilation of a circular loop. As an application in axion 

■ models, we discuss the relative loop and long string radiation contribution 

to a cosmological bound on the axion mass. A parallel linearized formalism 
for gravitational radiation is also provided. 



1 Introduction 

Global strings are associated with the breaking of a global symmetry that leaves 
a degenerate vacuum manifold which is not simply-connected. Global strings arise 
naturally in a variety of physical contexts; they are generic in axion models and 
seemingly also in superstring theories, and they can appear in GUT models, for 
example, through the breaking of a global B—L symmetry. For appropriate mass 
scales, global strings can take on many of the cosmological roles invoked for local 
gauge strings. Despite these similarities, the literature on global strings is sparse 
and their physics is inadequately understood and often even misunderstood. The 
chief distinction is that, unlike local strings which radiate gravitationally, oscillating 
global strings decay more rapidly into massless Goldstone bosons. It is the nature 
and characterization of this radiation which forms the subject of this paper. 

The popularity of the axion in particle physics and cosmology has provided one 
of the main motivations for studying global strings, so it is enlightening to provide 
a context for this discussion. The imposition of an extra global U{1)pq symmetry 
remains the most elegant solution to the strong CP-problem of QCD [PQ]. This 
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Peccei-Quinn symmetry is broken at some high energy scale and the correspond- 
ing pseudo-Goldstone boson 'd which appears is the axion [ Wil, Wei] . The axion field 
in these extended models becomes directly related to the previously arbitrary 9 
-parameter of QCD which is associated with any CP-violation. At the QCD-scale, 
"soft" instant on effects slightly tilt the symmetry breaking potential and the axion 
acquires a small mass ma- The minimum of the potential then becomes -i? = 
which corresponds to the CP-conserving value for 9. Note that neutron electric 
dipole measurements imply that 9 < 10~^°. 

In general, the couplings of the axion to ordinary matter are inversely propor- 
tional to the symmetry breaking scale fa- By making /a sufficiently large, the axion 
becomes effectively "invisible." A lower bound for /a is provided in astrophysical 
contexts by axion production because, being so weakly interacting, they can escape 
from the entire volume of a star and compete with surface photon or neutrino losses. 
A study of red giant cooling provides the limit /a > 10^ GeV [i?C], while supernova 
1987a implies f^>10^^ GeV [SN]. 

As one of the prime cold dark matter candidates, the cosmology of the axion has 
attracted considerable interest. Requiring that the density contribution of axions 
is below critical (that is, Qg, < 1) provides an upper bound on /a. If the universe 
passes normally through the Peccei-Quinn phase transition, then there are a number 
of possible sources of axions. Thermal axions produced by the phase transition are 
massless and are redshifted away by the expansion. Another source is the production 
of axions through oscillations about the minimum of the potential i!} — due to 
initial homogeneous misalignments of the axion above T > Aqcd- These zero- 
momentum axions were estimated to give the relative energy density, [v46',Pl^l^, 

"<.«o.9xio-=(^5A_).7 (1) 

where 9i is the initial average misalignment angle.* These axions were mistakenly 
thought to be the only significant contribution to the axion density, imposing the 
upper bound 

fa<10^^GeV. (2) 

This cosmological estimate ignored topological effects at the Peccei-Quinn 
phase transition [VE], that is, the inevitable production of global C/pg (l)-strings 
through the Kibble mechanism, f Davis [Dava,Davb] suggested that the radiative 
decay of the global string network provided a more significant contribution to the 

* Note that this bound can be weakened if there is substantial entropy produc- 
tion at some time between the QCD scale and nucleosynthesis, for example, by 
engineering the out-of-equilibrium decay of a massive particle species (such as ref. 
[Lythl]). 

t We should point out here that it is possible to avoid axion strings altogether 
through inflation, though such models may have unsatisfactory aspects. If /a is 
pushed above a GUT scale reheat temperature, we must by (2) anthropically arrange 
a small misalignment angle 9i in our region of the universe. The alternative of a 
low reheat temperature below 10-^^GeV assumes baryogenesis at the electroweak 



2 



axion density. This was estimated by assuming that direct radiation by long strings 
into the fundamental harmonic was the primary energy loss mechanism. Long 
string radiation was subsequently studied in more detail by Davis & Shellard [DSa] 
who demonstrated that strings radiate with a classical spectrum, predominantly in 
the second harmonic. Applying this spectrum with the previous assumptions and 
numerical results from network evolution, they deduced a cosmological bound 

/„<10i°G'eF, (3) 

which is two orders of magnitude stronger than the zero- momentum bound from (1). 
A comparison with the astrophysical constraint brought into question the viability 
of axion models in such scenarios, suggesting that inflation may be required to 
eliminate axion strings. 

There are, however, further contributions to the axion density. When the axion 
mass switches on at T = Aqcd, the 5"^ vacuum manifold is tilted, and the minima 
lie at the discrete values — 0, ±27r, ±4n, ... Thick domain walls form in regions 
where there are large variations in the axion field [Sika] . In particular, strings be- 
come attached to domain walls which begin to dominate the dynamics of the string 
network [ VE] . Reconnections between strings and domain walls facilitate the rapid 
decay and annihilation of the hybrid system into more axions, as demonstrated nu- 
merically [Sa,PRS,Sb]. These additional axions will serve to strengthen the bound 
(3). Paradoxically, axions produced by loops will tend to weaken the bound (3), 
which was derived assuming only long string radiation. If axions come predomi- 
nantly from small loops instead of long strings, they will be produced at higher 
frequencies and so lose more of their relativistic energy through redshifting [Sb]. 
We shall return briefly to this issue in Appendix A, leaving detailed analysis for a 
future publication. 

The orthodox classical understanding of string radiation assumed above has 
been challenged by Sikivie and coworkers who have raised a number of important 
issues [HarS,HagS,Sikb]. As we shall point out in §2, the global string appears to 
be a non-local object with logarithmic contributions to its energy density coming 
from the Goldstone boson field on all length scales above the string core width 
S. Implicit in radiation calculations, therefore, is the assumption that the global 
string maintains its integrity during its motion, and that these less localized compo- 
nents merely renormalize the string tension (provided that they are on length scales 
below the string curvature radius). Harari & Sikivie [HarS\ questioned this renor- 
malization and suggested that the string core during its motion somehow decouples 
from the surrounding "cloud" of Goldstone bosons. They proposed that the global 
string was in fact critically damped and lost all its oscillation energy in less than 
one period, implying that the radiation power spectrum was flat with contributions 
at all wavelengths above 5. These suggestions were supported by later numerical 

scale and may be ruled out in future by constraining the possible shape of the 
inflaton potential using microwave anisotropics, given that one knows the relative 
contribution of gravity waves [CKLL]. 
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simulations of perfectly circular loops which, by virtue of their symmetry, are in 
fact expected to behave in this manner [HagS\ . If verified, a fiat radiation spectrum 
would lead to a cosmological constraint, /a^lO^^GeV which is comparable to the 
zero-momentum bound (2) but at variance with original string bound (3). 

Here, we endeavour to settle this controversy beyond reasonable doubt; we 
make an extensive study employing sophisticated numerical techniques which we 
compare quantitatively with the analytic formalism developed to describe global 
strings. This formalism using two- index antisymmetric tensors [Wit, VVa,DSh] has 
been used to study radiation from loop trajectories [VVa,DQ] and periodic long 
strings [<S'aA;a]. We point out where our analysis differs from previous work and we 
introduce a simple backreaction model, based on the analytic methods, as a first 
step in describing the damping of small-scale structure on cosmic strings. Resolution 
of these questions is all the more pertinent because of ongoing and proposed axion 
dark matter searches. The motivation for these expensive experiments relies heavily 
on the theoretically predicted parameter range allowed by the cosmological bounds 
(2) and (3). 

Throughout this paper we employ a (-1 ) signature for the spacetime met- 
ric Qnui (H — ) for the induced metric on the string worldsheet 7a6, the coordinates for 
which are given by = X^{a, r), with the null coordinates, u = a — r, v = a + r. 

2 Analytic formalism 

The Goldstone model 

The essential features of global strings are exhibited in the simple U{1) Gold- 
stone model with action given by. 



s = J - ^A($$ - )2| 



(4) 

2\2 



where $ is a complex scalar field which has been split as $ = (/>e*'* into a massive 
(real) Higgs component (f) and a massless (real periodic) Goldstone boson -d. The 
corresponding Euler-Lagrange equations are given by 

d^d^^^+^X^^^-fi) = 0. (5) 

For a straight global string lying along the 2;-axis, we take the ansatz 

$(r,^) = 0(r)e-^ (6) 

where 6 is the azimuthal angle and n is the winding number. We take the usual 
boundary conditions, 0(0) = and — > /a as r — > cxo. Despite these conditions, 
the energy per unit length fi is mildly divergent (n=l). 



//(A) w //o + y --Q^ 2nrdr ^ + 2nf^lii \ — j 



(7) 
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where 5 ~ (V^/a)~^ is the string core width and no ~ /f is the core energy 
associated with the massive field cj) (that is, within r^S). The length-scale A is the 
renormalisation scale provided in general by the curvature radius of the string or the 
average inter-string separation. Typically, at the present day with /a ~ 10-^^GeV the 
logarithm in (7) is ln(A/5) ~ 70. This implies that there is much more energy in the 
Goldstone field than in the string core fiQ, a fact which has made the understanding 
of global strings intuitively hazardous. 

Dual representation: antisymmetric tensors 

The analytic treatment of global string dynamics is hampered by the topolog- 
ical coupling of the string to the Goldstone boson. However, we can exploit the 
well-known duality between a massless scalar field and a two-index antisymmetric 
tensor field B^^ to replace the Goldstone boson in (4) via the relation, 

cf>%^=lfae^.Xpd''B^P. (8) 

The transformation in the presence of vortices where — > has to be performed 
carefully [DSb] , (4) then becomes 

S = J cZ^x|(a^<^)2 - - iA(0^ - - 2nfaJ B^,da>^^ , (9) 

where the field strength tensor is Hf^"^ = d^^B"^ + d^Bf^" + d^B^^" and the area 
element rfcr^'^ is given in terms of the worldsheet X{a, r) swept out by the zeroes of 
the Higgs field (0 = 0), 

^^n,. = e''^daX''dbX''dadT . (10) 
Effectively, the string acts as a source term for with a current density given 

by 

ji^- = k j S{^)[x-X{a,T)\da'"' . (11) 

Integrating radially over the massive degrees of freedom in (9) yields the Kalb- 
Ramond action [ifi?], 

S = -fioJ dadr -^Jd^xH^- 2tt fa J B^^daf"" . (12) 

The first term is the familiar Nambu action for local strings, the second term is 
the antisymmetric tensor field strength for both external fields and the self-field of 
the string, and the last term is a contact interaction between the B'^'^ field and 
the string worldsheet. The coupling between the string and S^jy is analogous to 
the electromagnetic coupling of a charged particle, and is amenable to the same 
calculational techniques. 
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Varying the action (12) with respect to and B^^i, gives the equations of 
motion for the string and the antisymmetric tensor field, 

^JiQ{X'^ - X"^^) = 27TfaH^-^iX^X'p - X'^Xp) , (13) 
d'^daB^, = 27T fa J dadT{X^Xl - X;X,)5(4) {x - X{a, r)) , (14) 

where we have employed the conformal string gauge and the Lorentz antisymmetric 
tensor gauge, 

X^ + X'^ = 0, 

X.X' = 0, (15) 
d^B"" = . 

As for the electron we have to perform a renormalization because the antisymmetric 
tensor field comprises both the self-field and the radiation field of the string. The 
self-field is logarithmically divergent but it can be absorbed such that [LR,DQ] 

H{A){X'' - X"'') ^ f (16) 

where A is the renormalisation scale introduced in (7) and /'^ is the finite radia- 
tion backreaction force on the string. The exact form and effect of the radiation 
backreaction force term is the subject of another publication [BSb]. 

The renormalised equations of motion (16) for the string, assuming the effects 
of radiation backreaction to be small, that is /'^ ~ can be approximated by the 
Nambu equations of motion. Using the conformal gauge (15), these equations of 
motion have the general solution 

X^^t^T, X = I [a(w) + b(^;)] , (17) 

where u and v are null coordinates on the worldsheet and 

a'2 = 1, b'2 = 1 . (18) 

These equations have closed loop solutions and long (or infinite) string solutions. 
Fig. 1 illustrates two different long string solutions which are periodic and will be 
employed in our subsequent analysis. 

Radiation calculations using antisymmetric tensors 

Methods developed for calculating gravitational radiation from strings [Tur, VVb]^ 
[Bur] are readily applied to an antisymmetric tensor field radiation [FFa]. This 
work focussed on loop radiation, but the formalism has also been adapted for the 
study of long strings, that is, infinitely long periodic sources [Saka,Sakb]. Given a 
periodicity length L in both a and r, then along a 2;-directed string we can write 
X^ = aa + Pi{u) + P2{'v), where the functions Pi and P2 are also periodic with 
period L and assumed to be small relative to a. By analogy with the gravitational 
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Figure 1: Some specific periodic long string solutions: (a) helix, (b) pure 
sine wave 



radiation analysis [Saka], one can deduce that the average power radiated per unit 
length by the string is 



dP °° r"^-^ ~ 



"=1 \Km\<OlOJn 



where = 27rn/L, Km — 2nm/aL, = |k^|(cos(/?, sin*/?), |k^| = y/c^^"^^"^^ and 
J^'^ is the Fourier transform of the source distribution J'^'^ given by (11). Using 
the split into left- and right-movers for J^^'^ derived in Apppendix B one can deduce 
that 

-^T," H J iA\uf\vf-\u:vn, (20) 

n—1 |m|<c<n ^ 



dz 



m-\-n even 



where 



Jo L 

V^^= r -d^X^e-^'^-''-^^^ ^^^^ 
Jo L 



— (j^ni k , Km) 
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As described in the previous section, the antisymmetric tensor notation in principle 
aUows the calculation of the radiation power for particular solutions of the Nambu 
equations of motion. One important consequence of this formula is the fact that 
there is no contribution to the radiation power from the fundamental mode (n=l). 
(Note that the corresponding expressions for gravitational radiation power are given 
in Appendix C.) 

Linearized radiation power expression 

Unfortunately, it is still difficult in most cases to calculate the radiation power 
using (20) . However, this expression can be linearized by considering only first order 
terms in the modulus of radial perturbations [Hina\. One infinite periodic solution 
of (17) and (18) is 

X'' = (r, X^, aa + Px{u) + P2{v)) . (22) 

In linearizing, we shall ignore powers of l^"*"]^ or greater. The actual radiation 
expressions previously derived for gravitational radiation are not correct. Here we 
deduce the analogous expression for Goldstone boson radiation from a global string, 
while the corrected result for gravitational radiation can be found in Appendix C. 
The gauge constraints (15), in null worldsheet coordinates (tt, v), can be written 

as 

{duXY + {du^^f = i, {d.X^f + (a.X^)2 = i . (23) 
Rearranging and using a Taylor series expansion gives 

a.X^^ 1(1 -2(9^X^)2), a,X3^i(l-2(a,X^)2). (24) 

Hence to first order, X^ = (u + v)/2 so a — 1, Pi = and P2 = 0. We can write 
the linearized expression for Xr and Xl as 

= Xi(«),i«) X^{v) = {lv,Xi{v),lv). (25) 

These expressions can be inserted into (21) for and V^. Integrating U^, V^°, 
and by parts and linearizing the exponential then gives 



\27T{m + n) 27T{m + n) 

V'^=( rk^.V^, V^, ^k^.V^ 

\ 27r(n — m) 27r(n — m) 



(26) 



where 





' du 


/ 


~L 


/ 


' dv 


^0 


T 


lk^ 


■ n = 



i7r(m-\-n)u/ L 



(27) 
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Therefore, one can deduce expressions for \U\'^, I^Pj and \U.V\ in terms of 

the hnearized variables , and n, that is, 

= -|U^|^, = -|v^|^ 

t/*.F = 2(n.U^*)(n.V^) -U^*.V^, (28) 
U.V = 2(n.U^)(n.V^) - U^.V^ . 

We are now in a position to evaluate the (y?-integration as the only remaining depen- 
dence is in n. The following formulae for W,X,Y,Z (all two dimensional vectors) 
are easily obtained, 

/ d(/?(n.X)(n.Y) = ttX.Y (29) 
Jo 

/ dip (n.W)(n.X)(n.Y)(n.Z) = 

^0 



(W.X)(Y.Z) + (Z.W)(X.Y) + (X.Z)(W.Y) 

(30) 



Applying these formulae to the power expression (20) with (28) yields 

n=l |m|<7i 
m+n even 

If one of the components of X-*- is zero then the linearized power expression can be 
simplified further 

^ = ^Vn y lU^l^lV^r (32) 

n=l \m\<n 

m-\-n even 

The linearized power expression (31) will only give the leading order term if the 
average value of IX-*"! is sufficently small, note also that it is possible that there is 
no contribution at this order. Note also from (20) or (31) that left-moving modes 
propagating along a straight string will not radiate if right-moving modes are absent 
(and vice versa). 

Radiation power for specific solutions 

Long string configurations can be quantified by the wavelength, A (= L) and 
the amplitude, A. For solutions of the form (22), the important parameters are the 
ratio of the amplitude to wavelength, s = 27rA/X (denoted the relative amplitude), 
and the oscillation frequency, Q = 27t/X. In particular, we note that solutions 
for which e > 1 are not allowed by (18) in our linearized description. Of course, 
solutions with e > 1 are allowed for more general string configurations, but in the 
cases of high symmetry below they become highly degenerate with whole sections 
of string approaching luminal velocities (see the discussion in section 3). 
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Figure 2: Log-linear power spectrum of the helix solution at small (e = 
0.2) and large {e — 0.99) amplitude. 

(a) Helix solution 

A simple, symmetric solution of the form (22) can be constructed from equal 
and oppositely propagating sine waves, 



X = ( [ cos Qu + cos Qv] , — [ sin Qu + sinQv] , - 1 — e'^ {u + v)j 



(33) 



where < £ < 1 and e — > 1 in the relativistic limit. This corresponds to a helicoidal 
solution which oscillates between a static helix and a straight line. 

In this case we can use the full nonlinear analysis (20). The radiation power ex- 
pression contains many Bessel functions with argument k = \k^n^ — m? /{I — £^), 

f = E^" = E E P"."^' (3") 

n=l n=l i^K^yirp- 

m-\-n even 

where 

Pn,m ~ _ g 2 ~ \^^{{JlJu+l — JuJl+lf + {JuJl-1 — JlJu-lf — '^Jf Ju (^5) 
— \{Jl+lJu-l — Jl-lJu+l)^) + e^JiJu{2JlJu — Jl+lJu-1 — Jl-lJu+l) 
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and 

I — ^(m + n) , u = |(m — n) . (36) 

(This expression was previously obtained [Sakb] , but without the shortcut provided 
by the left- and right- moving split of Appendix B.) 

It can be easily shown that Pn,o = for all positive integer n. From this and 
the conditions \m\ < ny/l — and m + n even, one can deduce that the lowest 
even harmonic radiating is given by 



'^even = minjn : n > . .n an even integer} (37) , 

vl — £ 

and the lowest odd harmonic is given by 

nodd = minln : n > , n an odd integer} (38) 

V 1 — £^ 

For all £, neven is greater than two and, for e < Scrit = 0.94, nodd is three. Therefore 
for s < £crit) the third harmonic is the dominant contribution to the radiation power 
of the helix. In this case, one can deduce a linearized power expression by calculating 
the power in the third harmonic to first order in e, which gives 

dP 47r3/a^'° . X 

— ~ — . (39) 

dz 125L ^ ^ 

Applying the linearized power expression (31) to this solution gives zero, as it only 
takes into account terms of order and below. 

For s < Scriti one can see from fig. 2 that the power in the nth harmonic decays 
exponentially, with exponent proportional to n for large n (that is, oc e"""^, 
n » 1). Fig. 3 illustrates how a varies with e, showing that the exponential fall- 
off a becomes larger at small e. The proposal of Sikvie et al. for radiation from 
long strings requires that Pn oc n~^. This contention should also be contrasted 
with the result for the Burden-Tassie loop solutions {Pn oc n~^/^) and kinky loops 
{Pn OC n-2) [VVa]. 

The helix solution has a great deal of symmetry which probably causes the 
cancellation of the second harmonic. Configurations of this type would not be ex- 
pected to form in a realistic string network, because oppositely propagating modes 
will not have strong correlations. However, given that the helix is so weakly radiat- 
ing, it conceivably could be the endpoint of a radiation process, say for small-scale 
structure on a large loop. The stability of this solution to radiation damping and 
other perturbations is, therefore, an interesting issue. Given that we do not observe 
the cancellation of the second harmonic in more general solutions, we believe that 
the helix is actually an exceptional result and will be unstable to stronger radiative 
decay (a contention for which we provide numerical evidence later). Nevertheless, 
it is not unreasonable to suggest that some helix properties, particularly the expo- 
nential fall-off of the harmonics at large n, will be generic for all £ < 1 long string 
solutions. 
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Figure 3: a against e such that Pn oc e for large n. 
(b) Sine solution 

A solution, similar to the helix, but with perturbations in only one direction 

is 

X= — cos Qu + cos riv] , 0, — — - \_E{s, Qu) + E(^s, Qv)^ | , (40) 
where E{k,(f)) is the incomplete elliptic integral of the second kind, defined by 

E{k,(j))= r (41) 
Jo 

where < £ < 1 and £ — > 1 in the relativistic limit. 

The elliptic function can be written as a linearly increasing function plus a 
periodic function, so the solution can be seen to be in the form (22) and we can 
employ the power expressions (20) or (31). Unfortunately, the full nonlinear version 
(20) becomes too complicated, so we apply (31) to obtain 

— = . (42) 

dz 16L ^ ' 

This only yields the power up to order £^, for which only contribution comes from 
the second harmonic (n=2). Presumably, as for the helix, the power will become 
dominated by higher order terms at large £, but it seems reasonable to assume that 



12 



there also exists some £crit below which (42) is a good approximation. In any case, 
for small amplitude waves it is clear that the dominant contribution to the radiation 
power comes from the second harmonic. 

(c) Kink solution 

It is well known that (18) has solutions with points at which the string velocity x 
and tangent x' are discontinuous; these points, known as kinks, propagate along the 
string at the speed of light. Here wc give a simple solution for a periodic distribution 
of kinks on a long string, consisting of left and right moving perturbations, Xjj and 
Xy respectively. 



X= (Xc, + Xy,0,- 



(43) 



where 



Xu={ 



2e 



— u 

TT 

— ( -L+u 

TT 



< 


u < 


4^ 










< u 


<^ 






<L, 




< u 







(44) 
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2e 
— V 
TT 

2efl 

— -L — V 

n \2 
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V < 
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< V 


<^ 
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<L, 




< V 


4 





(45) 



where < s < ^ and, in this case, e — > is the relativistic limit. 

Once more the full non-linear radiation power expression (20) is intractable, so 
we use the linearized version (31) which yields 



dP 

dz 



2^4 



5L 



TT 



n=l 



I m I <n 
m-\-n even 



(n^ — m?)'^ 



1 — cos 



7r(m + n) 



1 — cos 



7r(m — n) 



(46) 

Again, this expression should be valid for s < £crit for some £crit) but even for 
sufHcently small s the radiation power is a logarithimically divergent infinite sum 
(that is, Pn oc for large n). This can be seen analytically from (46) since 
the sum is dominated by the m = n — 2 terms and numerically in fig. 4. The kink 
solution (43) is an infinite Fourier series over odd harmonics. However, in a physical 
context this divergence is cut off initially at the string width S and at increasingly 
larger scales as radiation backreaction rapidly 'rounds' the kink. This 'rounding' 
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Figure 4: Log-log power spectrum for the kink solution showing power 
law fall-off. 



will eliminate the more rapidly decaying high frequency contributions, ultimately 
leaving only the fundamental mode excitation. 

(d) Spectrum, of sinusoidal solutions 

The following solution corresponds to a superposition of a large, but finite 
number, of sinusoidal perturbations. 



(47) 



such that 



e 

mi 



AiNVLu I > ^ 



iNClv 



where Aj^ and Sjv satisfy 



(48) 



'-TV 



-A 



-AT, -DAT 



-B 



■TV 



An.A^-n = 0, Yl BTV.B„_Ar = Vn ^ . 

TV^O N^O 



(49) 
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Such a random superposition more realistically models the perturbations expected 
on a string network. Using the linearized radiation power expression (31), we obtain 



dz 2L 

n=l |?7i|<Ti 

m+rt even 



I -A. m-j-n I 



2|-D |2 

m — n 



2 2 



(50) 



I m-\-n . m — n I ~|~ I m-\-n ' m — n 
2 2 — 2 — 2 



This illustrates the generic £^ dependence of the radiation power even for the higher 
oscillation frequencies. In special cases, such as the helix, the contribution at this 
order can cancel. 

Simple backreaction model 

The string configurations discussed in the previous section are solutions of the 
Nambu equations of motion. In flat space they will continue to oscillate forever at 
constant amplitude. Even in an expanding universe, the amplitude will only decay 
slightly due to Hubble damping once A < H~^. This is clearly physically unrealistic, 
as it takes no account of the fact that such configurations radiate Goldstone bosons. 
Here, we will endeavour to incorporate this radiation backreaction in some special 
cases. 

In the previous section, we calculated the radiation power from several solu- 
tions, giving special consideration to the contribution from higher harmonics. This 
was only significant at large relative amplitude e or for 'kinky' solutions. Even in 
these situations, however, our expectation is that these higher harmonics will be 
rapidly suppressed as their amplitude decays or the kink is rounded. Radiation 
power is soon dominated by one particular mode. For backreaction in this simple 
case, the expressions for the radiation power include powers of e and one can deduce 
a phenomcno logical model for the decay of e. 

Given that the generic radiation power for any oscillation mode is proportional 
to we can model radiation damping by considering the following equations,* 

E = ajie^ , 
dP__dE_Pe^ (51) 
dz dt L ' 

where E is the energy per unit length and a is the solution-dependent coefficent of 
£^ in the power series expansion of 



^da. (52) 
dcr 



* Note we implicitly assume that the cut-off A for // in (51) is comparable to the 
curvature radius of the oscillating string. 
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Figure 5: Oscillation energy against time For A=10, B=100 (unbroken 
line) and B=1000 (broken line). 



For the sinusoidal solution a = |, whereas for the helix a = | (as well as /? = 0). 
In flat space, fx and L are constants so we can integrate the differential equation 
(51) to obtain 

£ = £o(^l + |^t) \ E = fx + afxsl(^l + ^?j . (53) 

where initially £(0) = Sq. Defining A = ji, B = a/iSQ and C = a/iL/PsQ, this 
becomes ^ 

E = A + b(^1 + ^'^ . (54) 

Initially, the decay is linear (t << C), while asymptotically at late times {t » C) 
the energy decays as a power law, t~^. Fig. 5 shows the oscillation energy against 
time for some reasonable values of A, B and C. A useful property of the model is 
the effective half-life of the oscillation energy, 

ti/2 = a/iL/Pel . (55) 

The half-life depends on the initial relative amplitude Sq, increasing rapidly as Sq 
decreases. 
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The power law decay (54) of the oscillation energy provides a useful quantitative 
test of this formalism (see section 3) , but it is a special case depending upon the left- 
and right-moving modes having equal amplitude. Strings in a realistic network will 
in general have a random superposition of unequal left- and right-moving modes. 
We can demonstrate heuristically that exponential decay will be generic for such 
modes. Suppose that and are the 'averaged' amplitudes of the left- and right- 
moving waves respectively. Motivated by (32) we might expect these modes to obey 
a cross-coupled version of (51) [Hinh], that is, 

^^M = _M£i, e^(0) = ^, 

,2 o 2 2 (56) 
«^^ = -^' -(«) = ^- 
The solution we find in this case when A ^ B is 



el = 



^2 _ j^2^-/3{B^-A^)t/aixL ' 

The asymptotic behaviour of these solutions is of interest. In particular, we find that 
for 5 > A, £^ — i> and ej^ ^ B"^ — exponentially, whereas for B < A, e\ ^ 
and e\ — B^. Whichever of the movers initially has the largest amplitude will 

approach a finite amplitude, while the other will decay to zero exponentially. The 
oscillation energy decay timescale again depends on the relative initial amplitudes. 

3 Numerical methods and results 

Algorithms and boundary conditions 

As we have seen there is a duality between the Kalb-Ramond action (12) and 
the Goldstone action (4) in the presence of vortices. This duality implies that the 
two actions should exhibit the same dynamics. In the previous section the dynamics 
of various long string configurations were examined within the framework of the 
Kalb-Ramond action. Here we directly compare the results of simulations with the 
same configurations in the Goldstone model. We tackle this problem numerically 
because, as we have noted, the Goldstone action does not lend itself to simple 
analytic calculations due to the topological coupling of the massless field. 

We have developed sophisticated numerical algorithms to allow the simulation 
of long string configurations within the framework of the Goldstone model. As in 
ref. [^c], we discretize space on a three-dimensional grid with dimensions A^i, A^2, -^3 
in the a;, y, z directions respectively, solving the rescaled (/a — > 1, A — > 2) Euler- 
Lagrange equation, 

^2^ Q-i^ 
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We employ a second-order leapfrog algorithm for the time derivative and second- or 
fourth-order finite difference approximations for the spatial derivatives (generally 
the latter). We denote the two algorithms as 2-2 and 2-4 respectively; they are 
discussed in greater detail elsewhere [BSa] (see also ref. [Sd\). 

The cylindrically symmetric straight string ansatz (6) for the field requires that 
we solve the boundary value problem, 



This was solved numerically using relaxation techniques. The finite boundary con- 
dition (f){r — > i?) = 1 introduces an error of order R~^, which we minimize using 

a large grid and spatial step-size. We can perturb the straight string ansatz using 
the solutions defined in §2 with relative amplitude e. We only impose perturbations 
in the x- and y- directions, thus introducing higher order errors, 0{e'^). However, 
these are negligible for small amplitude perturbations and so we effectively create 
highly accurate numerical solutions of the Nambu action. 

It is appropriate at this point to comment on the string self-field and the ra- 
diation field since they are not easily distinguished in the Goldstone formalism. 
There are, however, particular configurations in which the two can be separated. 
For example, when a perturbed string straightens the massless self-field becomes az- 
imuthal around the string as in (6) and the massive self-field can be easily calculated 
numerically from (59). In our simulations, this allows a straightforward subtraction 
at specified times to yield both the massless and massive radiation fields. Again, this 
procedure will be most accurate for small e, though it still yields useful information 
for £ ~ 1. 

Naturally the spatial discretization introduces boundaries into an otherwise 
unbounded problem. These artifical boundaries have an important effect on the 
numerical solution, especially at late times in problems involving radiation. For the 
long string configurations on the 2;-boundary we impose simple periodic, Dirichlet 
or Neumann boundary conditions. As in ref. [DSa], on the x- and j/-boundaries 
we employ absorbing boundary conditions, except here using a higher-order second 
approximation version [EM] [EH] [BJR] (refer to ref. [BSa]). These absorbing 
boundary conditions require the solution of 



at a; = and x = NiAx and a similar equation at j/ = and y = N2AX. These 
remarkable boundary conditions reflect a wave incident at 45° with only 3% of the 
original amplitude, that is, the reflected wave energy is reduced by a factor of 10""^. 
Fig. 6 illustrates the efficacy of these boundary conditions by contrasting the absence 
of reflected radiation from absorbing boundary conditions with large reflections 
from Neumann boundary conditions for the same t. These allow the simulations to 
be evolved for many grid light-crossing times without signiflcant influence on the 




(59) 



0(r ^ 0) = , 0(r 



00) = 1. 




(60) 
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Figure 6: Comparison of (a) Neumann and (b) second approximation 
boundary conditions for typical string radiation at late times 

interior string solution. Note that the grid face boundary conditions (60) must be 
appropriately augmented along the grid edges and at the corners [BSa]. 

Several runs were performed with Neumann boundary conditions on the x- 
and y-boundaries to demonstrate energy conservation. Energy was conserved to 
within 0.2% for both the 2-2 and 2-4 algorithms after several thousand timesteps. 
Unfortunately, only first-order spatial differences were used in calculating the energy 
density so no marked diiference was evident between the two algorithms. We should 
point out that the boundary conditions (60) only absorb massless radiation, but 
strings emit both massless and massive radiation. The massive radiation in (f) is not 
easily excited and it is highly suppressed at small amplitude. Furthermore, because 
of dispersion, it does not reach the boundary until late times. For reasonable values 
of e this massive radiation is barely detectable and has no discernible effect on string 
motion. 

We usually employed a grid with Ni = N2 = 200, = 50, and spatial and 
time step-sizes Ax = 0.7 and At = 0.3 (well within the Courant condition for 
these algorithms). Initial string configurations were demonstrated to be numeri- 
cally stable for over 3000 timesteps, both in the interior and at the boundaries. 
The algorithm was used to study periodic helix, sine and kink long string solutions 
and circular loop collapse. Throughout the simulations we monitored the amplitude 
of the perturbed oscillation e and the total energy per unit length E. We also ana- 
lyzed the massless radiation spectrum using the three-dimensional spectral analysis 
techniques described in Appendix D. 

Our numerical results will be seen to be at variance with those of Sikivie et al. 
[HarS,HagS,Sikb] for three primary reasons: 

(i) Firstly, in ref. [HagS\ they analyzed the total massless field spectrum, including 
the string self- field, which introduces non-radiating higher harmonics into their spec- 
tral analysis. Fig. 7 demonstrates this clear difference between the total massless 
field spectrum and the radiation spectrum analysed after subtration of the string 
self-field. Notice that the self-field completely dominates the radiation field at this 
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Figure 7: Spectrum for (a) the total massless field (b) the radiation field 



amplitude. Moreover, the fall-off in the higher modes for the total massless field 
spectrum is power-law {Pn oc n~^), whereas for the radiation spectrum it is much 
more rapid, being consistent with exponential. 

(ii) Secondly, in ref. [HagS\ they inadvertently concentrate on degenerate configu- 
rations with £ > 1 which do not correspond to the Nambu solutions presented in 
§2. These configurations are highly degenerate and collapse ultra-relativistically; 
whole regions along the string would achieve luminal velocities* but for radiation 
damping. An appropriate analogy is the perfectly circular loop solution which has 
a similar degeneracy and a corresponding pathological radiation divergence. Ad- 
ditional difficulties with these configurations are the inaccuracy of the ansatz (6) 
at large amplitude and the fact that our comparative analytic results only apply 

* These luminal string regions are known as 'lumps', and are distinct from 'cusps' 
where only a single string point reaches the speed of light. Unlike "lumps', cusps 
have a finite integrated radiation power [ VVb] . 
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Figure 8: Isosurfaces of the energy density for the evolving sine solution 
(£o = 0.6) shown after each period. Radiation damping causes a noticeable 
decrease in the oscillation amplitude. 

in the e < 1 regime. As it happens for small simulation scales, radiation damping 
for these "luminal" e > 1 configurations is so effective that higher harmonics are 
only marginally more excited. Nevertheless, the e > 1 results give the impression 
that the string is critically- or over-damped, whereas the £ < 1 configurations ac- 
tually oscillate — as expected — like decaying Nambu solutions. In any case, such 
degenerate configurations would not be expected in a realistic string network, if 
only because of the presence of small-scale structure, and they are not, therefore, a 
reliable indicator of how it will radiate. 

(iii) Thirdly, the reader is invited to compare the rather different spectral analysis 
methods of ref. [HagS\ with those detailed in Appendix D. In particular, notice 
that the standard definition of the one sided three dimensional power spectrum 
(D3) and the related wave number (D2) are not used in ref. [HagS\. This apparent 
bin mislabelling will have an important effect on the numerically estimated power 
spectrum. 

Massless radiation from the sine solution 

The evolution of the sine perturbation (40) was simulated for a number of 
different initial amplitudes, {sq — 0.2, 0.3, .., 1.2) in order to obtain a quantative 
understanding of the evolution. As expected we find that the string oscillates losing 
amplitude by emitting radiation. Fig. 8 shows the decay of the maximum amplitude 
e. Emitted radiation is illustrated in fig. 9 in a plane transverse to the string at the 
point of maximum amplitude. The radiation pattern is dominated by quadrupole 
lobes of the n = 2 harmonic. 

We have analyzed the power spectra for eq = 0.4, 0.6, 0.9; these are illustrated 
in fig. 10. The subtraction of the string self-field from the much smaller radiation 
field has to be performed precisely. Nevertheless, we found that some minor resid- 
uals could not be eliminated, notably in the third and fourth harmonics, because 
of the numerical inaccuracies inherent in this procedure and the nature of the ra- 
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Figure 9: Radiation quadrupole patterns in the plane perpendicular to 
the string for eq — 0.6 after (a) t — 3X and (b) t — 7A. The difference the 
the radiation field form zero determines the colour saturation. 



Figure 10: Radiation power spectrum for the sinusoidal solution on a 
log-linear scale for (a) Sq = 0.4, (b) Eq = 0.6 and (c) Eq = 0.9. 

diation pattern near the string.* Configurations producing the strongest radiation 
provide the cleanest spectra {eq = 0.9 in fig. 10c). We can only be confident that 

* These appear to be residuals because they remain when the oscillation energy 
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Figure 11: The contribution of the second fundamental mode to the ra- 
diation spectrum for eq = 0.6. 

these results provide an upper bound for the radiation in higher harmonics (n > 1), 
but we can still draw definite conclusions. First, it is clear that the lowest radi- 
ating harmonic {n—2) provides the dominant radiation component, as predicted 
analytically. Secondly, although some higher harmonics appear to be present, they 
are strongly suppressed by radiation into the lowest harmonic — even for £o~0.9 
they are not significantly excited. This is an important observation which confirms 
the assumption underlying our simple backreaction model in section 2. Finally, 
the decay in the higher harmonics is consistent, within the uncertainties, with an 
exponential fall-off. Without doubt it is completely inconsistent with a flat power 
law spectrum, Pn oc n~^. An fall-off would predict an amplitude for the n=10 
harmonic relative to n—2 which is nine orders of magnitude larger than that ob- 
served in fig. 10c. The time evolution of the contribution to the radiation in the 
second harmonic is illustrated in fig. 11. The energy increases until the radiation 
propagates out to the edge of the box and begins to fall as the absorbing boundaries 
take effect. 

and the n=2 radiation pattern have largely died away. We believe that they are 
due to an imperfect subtraction because the strings is not exactly straight at the 
discretized sampling timestep. The effect is also exacerbated by some distortion 
due to Lorentz contraction and the fact that the radiation pattern is only purely 
quadrupole asymptotically, not in the immediate vicinity of the oscillating string. 
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A 


B 


C 


0.3 


28.18 


0.64 


769.02 


0.4 


28.14 


1.08 


386.45 


0.5 


28.08 


1.68 


238.68 


0.6 


28.02 


2.34 


165.04 


0.7 


27.98 


3.12 


123.91 


0.8 


27.90 


3.96 


99.37 


0.9 


27.84 


4.88 


84.20 


1.0 


27.78 


5.84 


74.82 


1.1 


27.66 


8.00 


65.81 


1.2 


27.50 


10.42 


64.73 



Table 1: Best fit values for the backreaction parameters A^ B and C. 



Figure 12: Energy against time: Numerical results and best fit curves 
for (a) eo = 0.5, (b) sq = 0.8 and (c) eo = l-l- 
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Evidence for quantitative agreement with the analytic radiation calculations 
comes from attempts to fit the time history of the overall energy density to the 
simple backreaction model (53). The total energy initially remains constant until the 
radiation propagates out to the boundary. At this stage it begins to fall with a time- 
delayed profile similar to that of our simple model. The least squares estimators 
for A, B and C in (54) are calculated using a global minimisation routine and are 
displayed in Table 1 for a variety of initial amplitudes. 

We find that A is constant to within a few percent and the calculated values for 
B and C are in good agreement for £o~l-l- The numerical results and the best fit 
curves are shown in fig. 12 for = 0.5, 0.8, 1.1. A more significant test comes from 
determining the dependence of B and C on ^o- This is shown in fig. 13 in a log- log 
plot, the points are well-fitted by a straight line in the region 0.3<£o^0.8. If we 
assume that A = Aq, B = BqSq^ and C = CqEq'^ then we calculate least squares 
estimators, Aq = 28.0±0.4, log^o = 1.7±0.2, log Co = 4.2±0.4, m = 2.0±0.1 and 
77,2 = — 2.0±0.1. These exponents are in excellent agreement with the predictions 
of the backreaction model (54). From these values we can obtain the parameters in 
(51), a — 0.2 ± 0.05 and (3 — 3.5 ± 1.75. Given the uncertainties, these are again 
in satisfactory agreement with the actual analytic values, a = 0.25 and /? = 1.98 
from (42). We can also calculate the value of the unrenormalized string tension by 
summing the massive field contributions to the energy of a straight string solution, 
yielding //q = 4.2 ± 0.2. This allows us to estimate log(A/5) = 3.8 ± 0.2. This 
indicates that the simulated global strings have approximately six times more energy 
in the Goldstone field than in the massive string core. These results, therefore, test 
and confirm the renormalization procedure which underlies the previous analytic 
predictions. 

We have also numerically studied the dynamics and radiation spectra of sinu- 
soidal solutions for eq > 1 where the analytic results are invalid. These solutions 
are highly degenerate and were expected to produce a large amount of radiation in 
higher radiation modes, however, this was not found to be the case. Fig. 14 shows 
the radiation spectrum for = 2.0, showing that the string predominantly radi- 
ates in the second harmonic even at these large amplitudes. This result indicates 
that radiation backreaction damps the excitation of higher harmonics much more 
strongly than expected. On these small simulation scales, massive radiation is also 
excited at large amplitude (see below) and this may contribute to the radiation 
damping. 

Finally, we should comment on the related helix solution which is predicted 
to radiate extremely weakly (39). Interestingly, we observe that the helix appears 
to be dynamically unstable and radiates in a similar manner to the sine solution 
(mainly in the n = 2 harmonic). This may be because our initial ansatz is not a 
sufficiently perfect representation of the true solution. However, we believe that the 
non-radiating property of the helix is a finely balanced result of its symmetry. Any 
asymmetry introduced by radiation damping or otherwise will tend to destabilize 
the solution towards others which generically radiate more strongly. This conjecture 
remains to be justified analytically. 
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Figure 13: Loglog plots of (a) B agianst Eq and (b) C against £o- 

Massless radiation from the kink solution 

We have also studied the evolution of the kink solution (43) for a number of 
different initial perturbation amplitudes. In section 2 we proposed that the kink 
solution should radiate with an spectrum, truncated at some time-dependent 
'rounding' scale. The effectiveness of this 'rounding' can be seen in fig. 15 which 
shows an initial kink perturbation evolving into a sinusoidal perturbation. The 
decay of high harmonics in the kink spectrum takes place very rapidly. Fig. 16 shows 
the radiation spectrum at early times (a), contrasted with the spectrum from late 
times (b) . Note that the kink solution initially has more power in higher frequency 
radiation modes with a slow power law fall-off. The later spectrum (fig. 16b), 
however, is indistinguishable from a typical sinusoidal spectrum. These results 
certainly support our intuitive understanding of kink 'rounding'. 

Massive radiation 

We have already noted that global strings can radiate massive particles, as 
well as massless particles. Depending on the final decay products, massive radiation 
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Figure 14: The radiation spectrum for a large amplitude sine solution 
with So = 2.0. Notice the power in the higher harmonics which shows up 
in the log-linear spectrum. 

from topological defects is of considerable interest because it has been proposed as a 
potential source of baryon asymmetry. Massive modes will be suppressed except in 
regions of high curvature, such as cusps and sharp kinks. In our simulations because 
the curvature radius is not much larger than the string width, massive radiation 
could be excited for large amplitude sinusoidal solutions, that is, £o~0.6. Massive 
radiation produced by global strings exhibits very different behaviour to massless 
radiation. Fig. 17 illustrates a circular pattern reminiscent of dipole radiation, quite 
unlike the quadrupole patterns of fig. 9. As the string oscillates, massive radiation 
is emitted most strongly in the direction of string motion. However, note that for 
amplitudes below eo < 0.9 the massive radiation signal was relatively very weak. 
Massive radiation will be discussed in greater detail elsewhere [BSa] . 

Loop solutions 

The radiation spectrum from oscillating loops was not easily accessible to nu- 
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Figure 15: Energy density isosurfaces for the kink solution [eo = 0.6) 
shown after each period of its evolution. Note the effectiveness of 'round- 
ing' due to radiation damping. 

merical simulation because on these small scales radiation damping leads to their 
rapid demise. There are further numerical difficulties setting up appropriate initial 
configurations and distinguishing the radiation field from the loop self-field. Nev- 
ertheless, we illustrate the collapse of a perfectly circular loop in fig. 18. The loop 
begins to shrink rapidly under its own tension, and annihilation occurs when the 
opposite segments collide. At this stage massive radiation is also produced copi- 
ously. The massless radiation pattern after collapse is shown in fig. 19. The fact 
that the radiation is maximal in a direction transverse to the loop plane is consis- 
tent with analytic calculations for the circular solution. The predicted divergence in 
the the integrated power is, of course, removed by radiation backreaction and finite 
width effects. We note that this very degenerate loop configuratoin is expected to 
produce an atypical power spectrum, so we cannot infer general conclusions for 
loop radiation from this special case (as in ref . [HagS\ ) . 

Our expectation on cosmological scales for non-intersecting loop solutions is 
that the antisymmetric tensor formalism of rcf. [FVa] will provide an accurate 
picture of loop radiation, that is generically oc n~" with a > 3/2 for n » 1. In 
the final stages of loop decay when its radius approaches the core width, the demise 
will proceed more rapidly through annihilation as in fig. 18. The resulting small 
burst of massive radiation may have significant cosmological implications. 

4 Discussion and conclusions 

The conclusions of this work are unambiguous. We have made careful compar- 
isons between high precision numerical simulations of the Goldstone model (4) with 
analytic calculations using the antisymmetric tensor formalism of the Kalb-Ramond 
action (12). The description of vortex- line dynamics with both approaches is found 
to be in close quantitative agreement. This has a number of important implica- 
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Figure 16: The radiation spectra for the kink solution after (a) t = 3A 
and (b) t = 15 A. Compare the latter with the sinusoidal spectrum given 
in fig. 7b. 

tions, not least because the validity and usefulness of this duality had previously 
been questioned. 

Firstly, this work confirms the procedure employed for passing between the 
Goldstone boson field and the antisymmetric tensor field B^'^ in the presence of 
vortices. Note that on the relatively small scales we consider in these simulations, 
the duality is not straightforward because the relation (8) is complicated by the 
behaviour of the massive field (f) in the string core. Provided perturbation wave- 
lengths are large compared to the string core width 5, we can integrate out the 
massive degrees of freedom to obtain the Nambu term in (12). Of course, for very 
high energy perturbations it proved possible to excite these internal massive degrees 
of freedom. 

Secondly, the Nambu-like behaviour predicted for global strings in (16) depends 
upon the renormalization of the string self-field. (This is distinct from renormal- 
ization about a point source because it involves cut-offs at both small and large 
scales.) In our numerical simulations, we typically had an energy in the string self- 
field (the Goldstone mode) which exceeded that in the core (the massive mode) by 
a factor of about six. Nevertheless, but for radiative damping, the evolution of ini- 
tially perturbed solutions was observed to closely correspond to oscillatory Nambu 
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Figure 17: Massive radiation contours for a sine solution with £o = 0.9 
after (a) t — 5\ and (b) t — lOA. Note that this is a weak signal with 
contours lying in the range, — 0.002 <1 — 0< 0.002. 



Figure 18: Energy density isosurfaces for a collapsing circular loop shown | 
at intervals of 20 timesteps. 

trajectories. In no sense is this self-field a diffuse 'cloud of energy' somehow distinct 
from the global string core. On the contrary, it is tightly bound to the string — the 
nearer the core, the more tightly. On scales below the string curvature radius, these 
massless modes merely renormalize (and dominate) the string tension and energy 
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Figure 19: Massless radiation contours in the aftermath of circular loop 
annihilation. Illustrated is a slice transvere to the plane in which the loop 
collapsed; very little radiation is emitted in the loop plane. 

density. 

Thirdly, as predicted by the Kalb-Ramond action, global strings radiate with a 
classical spectrum in harmonics directly related to the wavelengths of perturbations 
on the string. For example, for the sine solution we confirmed the absence of 
the fundamental mode, demonstrating that radiation was primarily in the second 
harmonic. This discrete spectrum appeared to be consistent with exponential fall-off 
at high harmonics for the long string configurations we studied, in complete contrast 
to the prediction of power law fall-off, Pn oc n~^, suggested in [HarS,HagS,Sikb]. 
Sharp kink configurations initially showed additional power for large n, but this was 
rapidly eliminated by rounding due to radiation damping, and the spectrum soon 
conformed to that for simple sine solutions. 

Finally, we have numerically confirmed the quantitative rate of radiative back- 
reaction predicted for particular long string configurations. Both the amplitude 
and decay rate of the oscillations on the strings were found to be in good agreement 
with our analytic results. Moreover, we observed that radiation damping was very 
effective at eliminating higher harmonics n > 2, markedly increasing the ampli- 
tude range over which our simple backreaction model gave satisfactory results. We 
are currently applying this work to more general string configurations because it 
has important implications for the understanding of string network evolution and, 
specifically, the build-up and damping of small-scale structure on strings. Modifi- 
cations to (16) and direct comparisons with field theory simulations are the subject 
of a forthcoming publication [BSb]. 

It is interesting to place these results in a cosmological context because this 
formed our initial motivation. For heavy cosmological strings, the logarithm in /i(A) 
is much larger (ln(A/5) ~ 70-100) than for strings in our simulations. Although 
the emitted power remains the same, the overall decay rate is slower because the 
radiation damping term in (16) is relatively weaker. Conversely, however, these 
strings can be expected to produce more radiation in higher harmonics because the 
string trajectory will more nearly approach the strongly radiating configurations 
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allowed by the Nambu equations, notably cusps. 

We believe that these results leave little doubt that the low-energy effective 
Kalb-Ramond action provides an excellent quantitative description of global string 
dynamics and radiation. Of course, as we have observed with massive excitations, 
the antisymmetric tensor description needs to be supplemented with the full non- 
linear field theory in very high curvature regions, such as at points of string inter- 
section or cusps. However, these regions typically have length-scales corresponding 
to the string core width S, so in a macroscopic context they will have a very minor 
effect on overall predictions of massless radiation. In cosmology these results imply 
that the antisymmetric tensor formalism should provide an accurate picture of ax- 
ion production by an axion string network as described in ref. [VVa,DSa,Sakb], but 
contrary to the assertions of refs. [HarS,HagS,Sikb]. Whether or not this has been 
achieved with sufficient precision at this stage remains an open question which we 
briefly take up in Appendix A. However, axion strings are by no means the only 
application of this elegant and tractable formalism; global strings also appear in 
low-temperature and other contexts. Indeed, refs. [DSc,GKPB] demonstrated the 
close relationship between global strings in the Goldstone model (4) and vortex-lines 
in a 'relativistic' superfluid. In ref. [Car] this has been generalized to compressible 
superfluids and perfect fluids, so it may be that the methods described here — both 
analytic and numerical — can be usefully applied and tested in a laboratory setting. 
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Appendix A: Axion string bounds revisited 

The astrophysical and cosmological constraints on the axion were introduced 
in section 1. Here we consider the axion string bound in the light of our previous 
discussions, specifically noting the relative axion contributions radiated by small 
loops and long strings. Of course, we assume the standard axion scenario with the 
universe passing through the Peccei-Quinn phase transition after reheating. 

Cosmic string evolution in an expanding universe has been extensively studied 
numerically [AT,BB,AShe]. A long string network loses energy through loop pro- 
duction and is observed to evolve toward a 'scaling' regime in which its large-scale 
properties remain constant relative to the horizon. These robust large-scale results 
were obtained for local gauge strings but they can also be expected to apply to 
global strings. The overall long-string energy density, however, must be modified 
with a logarithmic time dependence arising through the cutoff (A ~ t) in the global 
string energy density fx, that is, 

Poo = ^«^^ln(tA), (^1) 
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where ^ ~ 14 [BB,AShe]. The network correlation length is then given approxi- 
mately by ^ = C~^^H. 

In numerical simulations, evolution on small scales does not actually 'scale', 
though this appears to have little effect on the overall density (Al) and correlation 
length ^. In any case, in a realistic setting, this small-scale structure will be limited 
by radiation backreaction, damping out perturbations below a constant relative 
scale, Lr = Kt. The long-string backreaction scale should provide an lower limit 
to the typical 'scaling' size of loops created by the network, £ = at. 

Axions will be radiated by loops and long strings until the axion mass 'switches 
on' at the time tw = ^qcd when domain walls form, causing the rapid demise of 
the network. We can provide a greatly simplified estimate of the resulting number 
density of axions Ua by noting that the dominant contribution comes from the least 
relativistic axions created just before tqcD [Davb] (refer to ref. [VS\). To maintain 
scaling in each Hubble time, the long-string network must lose most of its energy 
to loops which, in turn, rapidly decay into axions of roughly the same wavelength. 
Just prior to t^, then, loops will have made the following contribution to the axion 
energy density, 

Pa^~^^ln(W^). {A2) 

These axions will typically have a frequency, u> = An /at, because loops of length 
£ ^ at radiate primarily in their fundamental mode 2/£. For massless axions we 
have Pa = iT'g.u!, so the number density of axions remaining at is approximately 

n,i-^Ht^/S). {A3) 

Clearly, a determination of the loop-size parameter a is essential for a quantitative 
estimate of the axion density today. 

Loop creation will occur on scales limited by the long-string radiation back- 
reaction scale, that is, a'^K. Recalling the damping 'half-life' ti/2 = a'fiL/(3eQ 
defined in section 2, we can determine the minimum length-scale L k, nt which will 
have lost most of its energy in one Hubble time: 

«~^[ln(t/5)]-^ , {AA) 

where we have taken jS = and a' = 0.25 for the sine solution (47), while 

also assuming Eq ~ 1. This yields k ~ 0.02 at tqcD, but large uncertainties must 
be admitted. It seems reasonable, therefore, along with previous authors to take 
a ~ 0.1 [VVa,VS\. 

By equal matter-radiation teq, the now- massive axions in (A3) are non-relat- 
ivistic, so the axion density is given simply by pa ~ Wana, where Ug, has been 
appropriately diluted by the expansion a~^. A comparison at teq between pa and 
the critical density pc yields the cosmological constraints [ VS\ , 

ma>10-SV, /a<10i^GeV, (^5) 
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where we have assumed a ~ 0.1. 

A rough estimate of the axion contribution from long string radiation is also 
straightforward. According to (42), the power radiated per unit length is given by 

— = — — , (A6) 

where we assume perturbations are peaked at some length-scale 7t (k < 7 < C"^''^); 
for which the relative amplitude £ 1. We now multiply by t and the total length 
of string (/t^ to obtain the axion energy density radiated in one Hubble time prior 

to tyf, 

Pa.-^. iA7) 

Given radiation in the second harmonic co = Aiv/'^t, we obtain the long string axion 
number density 

naoo-^^^. (^8) 

64tw 

Note that the axion (or graviton etc.) number radiated by long strings appears to be 
independent of the actual perturbation length-scale which dominates the spectrum. 

A comparison between contributions from loops (A3) and long strings (A8), 
yields the relative density ratio 

- ln{t^/S) ^ — . (^9) 



If loops are created above the backreaction scale a^K, then the axion density 
produced by long strings is probably subdominant. However, if a ~ k these densities 
may be comparable. As mentioned in section 1, there are also other sources such 
as axion emission during domain wall decay. Again, approximate estimates suggest 
that this contribution cannot be neglected in determining a final cosmological bound 
on the axion [Lyta]. 

It remains to comment on the axion string constraints presented in earlier 
work. It is apparent that the bound (5) is weaker than the original string bound 
(3), /a^lO-^°GeV, of refs. [Davb,DSa]. The difference lies not in the understanding 
of global string radiation — which agrees with this analysis — but, rather, in whether 
the dominant network energy loss comes from small loop creation or long string 
radiation. The assumption [Dava] that it comes from long strings radiating at the 
correlation length ^ ~ t seems to give an order-of-magnitude overestimate (compa- 
rable to taking the loop result (A3) with a 1). The apparent severity of this old 
long string constraint is illustrated in fig. 20. 

The difference between (A5) and the string bounds of refs. [HarS ,HagS ,Sikh] 
are more fundamental because they are based on the assumption that global strings 
radiate with a flat power spectrum, P„ oc . This is contrary to our results 
which show that power is generally dominated by the lowest allowed harmonic. 
The remaining literature on cosmological axion constraints has been flawed either 
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Figure 20: Schematic of the approximate astrophysical and cosmological 
bounds on the axion. The viabihty of the axion in the standard scenario 
clearly depends on an improved understanding of axion string radiation. 
The conclusions of this paper indicate that the true constraint will lie 
somewhere near the string loop bound (5), leaving a a somewhat smaller — 
but still substantial — parameter window. 

by prevarication on this issue or by ignoring topological effects altogether. Authors 
have preferred to reserve judgment on the underlying global string spectrum, while 
discussing the consequences of both alternatives (see, for example, refs. [Lyta,KT]). 

This paper only constitutes a first step towards a definitive axion string con- 
straint. While it supports the assumptions on which (A5) is based, we make no 
suggestion that axion detection experiments should either be launched or cancelled 
on the basis of the uncertain parameter window shown in fig. 20. However, we 
do hope to remove these remaining uncertainties with large-scale string network 
simulations which incorporate radiation backreaction. 

Appendix B: Power expression split with left- and right-movers 

One can split arbitrary solutions of the Nambu equations into left- and right- 
moving waves by employing null coordinates on the string worldsheet u = a — r and 
V = a — T (ie. = + X^). This was first used by Burden [Bur] to separate 
the Fourier transform of the energy-momentum tensor into left- and right-moving 
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components when calculating the gravitational radiation power from a closed loop 
(see a more extensive discussion in [AShe2]). This split has also been performed for 
a long periodic string [Hina]; unfortunately, however, there was a slight oversight 
which has important consequences for the allowed spectrum of radiation. 

Here we outline the left- and right-splitting of the Fourier transform of the 
source distribution J^'^ for a z-directed long global string, periodic in a and r with 
period L and — aa (conditions satisfied by the helix solution (33)). We then 
discuss the generalization to the case where X^ = aa + Pi{u) + P2{v) such that 
Pi and P2 are arbitrary (small) functions with period L. In the special case of 
= aa, we have 

p POO /"OO 

J'^^(x^,^,t)=^ / da driXf'X"' - X'f'X'') 

^ J —00 J —00 

x5'^{^^ -X)5{z-aa)5{t-T), (Bl) 

J''''{k^,Km,OJn) £ dt j^"^ dz J rf^x^ J'^^(x^, t)e-^^-^ , 

for n 7^ and = (k-*-. Km, <^n) where cvn = 27rn/L , Km = 27rm/aL. Using these 
definitions one can deduce that 

>^(k^,«^,a;n) = -^ ^da ^ dT{X^X'^ - X'>^X^)e-'^'' , 
2ai^ Jo Jo 

= / da j drl. 
Jo Jo 



We now change the variables of integration from (cr, r) to the null coordinates {u, v), 
in which / can be re-expressed as 



X^X'" - X'^X- ^2{dyXlduX>'j, - duX'^^d^Xl^ . 



(B3) 



The region of integration is more complicated. Fig. 21 shows how the integration 
over the region Q < a < L and < r < L (that is, -Ri U -R2 U -R3 U -R4) can 
be transformed into an integration over the region < w < L and < v < 2L 
(i?i U i?2 U i?3 U i?4). This can be further simplified by noting that 

I dudvl = e'^'^+''^'' [ dudvl (B4) 

Jr4UR'2 JriUR'^ 



Therefore we can write the Fourier transform of the source distribution function in 
terms of left- and right-movers as 



J>^(k^, Km,U;n) = ^ [U'^V - V^U"] (B5) 
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Figure 21: Equivalent regions of integration for (cr, r) and {u,v). 



where 




— {(^m k , Km) • 



where m + n must be even. This condition was missed in the previous derivation 
and has a important effect on the allowed spectrum of radiation. In particular it is 
responsible for there being no contribution from the fundamental harmonic [n—l). 

The generalisation of this splitting to a solution with = aa + Pi{u) + P2{v) 
can be achieved, assuming the Pi and P2 are sufficently small by allowing the lines 
bounding the regions i?i,i?2,-R3 and R4 to become oscillatory. The periodicity of 
fig. 21 is maintained by the periodicity of Pi and P2- 



Appendix C: Formulae for gravitational radiation power 

This appendix gives gravitational radiation power expressions for local strings, 
analogous to the those for Goldstone bosons from global strings. The source term 
for the gravitational radiation is the energy-momentum tensor 

T^"" = 1^0 [ (X'^X" - X'^'X"')S^'^^ [x - X{a, r)] dadr , (CI) 
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and the radiation power is given by 

n=i m st \Km\<0('jJ 



J2 ^ d^(^T'^"*(a;n,k^,Av™)VK,k^,«r 



(C2) 



"m) I 



On performing the spht into left- and right- movers, as in appendix B, one can 
deduce that 



dP _ 327rG'/i^ 
dz 



^^f^n Yl r dif{\U\^\Vf + \U\V\^-\U-Vf}, (C3) 

n=l \7n\<otn 



m+Ti even 



and the corresponding hnearized power expression is given by 

dP _ GAnGul °° 



dz 



5^ {|U^|2|V^|2-|U^*.V^|2 + |U^.V^|2}. (C4) 



71=1 \m\<n 
m+rt even 



Note the difference between the previous expression (C3) and the hnearized version 
(C4) (a correction for ref. [Hina]). Applying (C4) to the helix, sine, kink and sine 
spectrum solutions discussed in section 2 yields the following 



Helix: 

Sine: 

Kink: 



dp 




dz 


L 


dP _ 


n'^Gfils^ 


dz 


2L ' 


dP 




dz 





(C5) 
(C6) 



■n= 1 I I < n 

m-\-n even 



X 



1 — cos 



7r(m + n) 



1 — cos 



7r(m — n) 



(C7) 



dP _ An^G{ile^ 
dz 



Sine spectrum: 



E E 

n=l |m|<n 

m-\-7i even 



\\ |2|-R |2 

x\. m-j-n Xj m — n 
2 2 



B 

2 2 



— I A* R 



(C8) 



This behaviour is very similar to antisymmetric tensor radiation. The notable 
exception is the helix solution which radiates gravitationally in the second harmonic, 
underlining the incongruous nature of the Goldstone boson result. Note also that the 
linearized power expressions above are only valid for e < £crit for some £crit = C'(l). 
A simple backreaction model for gravitational radiation will take the same form as 
(54). 
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Appendix D: Three-dimensional spectral analysis 

Standard methods for power spectrum analysis [W^,iJ,A^ are applied to the ra- 
diation field on a three-dimensional grid. This involves calculation of the three- 
dimensional discrete Fourier transform using an FFT algorithm. The three-dimen- 
sional Fourier transform, Fp^^p^^p^^ of some discrete function -Fa 1,02, as is given by 



-^/3i,/32,/33 - 

where 



^ ATi-l 7V2-I 7Vs-l / s 

* ai — 1 a2 — l a3=l 



Fai,a2,a3 - F{aiAx,a2Ay,a3Az) , 
F^uP.,P3 j^, ) , (D2) 

< ai, Pi < Ni - 1 for i=l,2,3 . 
Given the above, the one-sided power spectrum is defined as [PFTV] 



Par,a2a3 (^^^^^^^^2 



1-^01,02,0:3 I + l-^ATi— Oi,iV2— 02,-/V3 — 0:3 I 



(D3) 



for 1 < CKj < — 1) i=l,2,3. In the special case when one or more of the ctj is 
either zero or then adjustments must be made; in particular. 



2 

2 



2 ~ 2 

PjVi JV2 £3 = , \Fnj^ N2 jVg 



2 

(D4) 



2 ' 2 ' 2 iViiVoiVc!) 2 ' 2 ' 2 



and similar expressions are used for ^0,02,03 ; -^0,0,03 etc. Once we have the three- 
dimensional one-sided power spectrum, we calculate the modes to which the power 
spectrum corresponds. This is done by comparing the wavelength of the string with 
the wave number of each point in the power spectrum and putting the contribution 
into an appropriate bin. 

There is a problem associated with this approach to power spectrum analysis. 
The discrete Fourier transform assumes that the initial data is periodic, which is 
not necessarily the case for the radiation field in our simulations. The resulting 
discontinuities at the boundaries can cause spurious harmonics to appear in the 
spectrum. We overcome this by multiplying the data by a window function, W{x). 
Window functions should be periodic and fall-off quickly toward zero at the edges. 
To apply the window function we calculate the discrete Fourier transform of G{x) = 
W{x)F{x). Using the convolution theorem and the properties of W{x), we can 
recover the Fourier transform of F{x) as 

f3i,(32,(33 - > K^O) 
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where 




N1-IN2-IN3-I X 2 



(D6) 



The window function VFq,^ ^0,2^^3 is decoupled in each direction, that is, Wa-i^^a2,a-i = 
Wctj^Wct^Woi^, where each of the W^^ is a one-dimensional window function similar 
to those described in ref. [H]. After extensive trials, it appeared that the problem 
was best suited to a Gaussian window function, that is, 



where i — 1, 2, 3 and (3 is comparable to the inverse of the simulation box-length. 
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